GLOBAL EXISTENCE FOR THE VLASOV-POISSON 
SYSTEM IN BOUNDED DOMAINS 

HYUNG JU HWANG AND JUAN J. L. VELAZQUEZ 

Abstract. In this paper we prove global existence for solutions of 
the Vlasov-Poisson system in convex bounded domains with specular 
boundary conditions and with a prescribed outward electrical field at 
the boundary. 



1. Introduction 

In this paper we study global solutions for Vlasov-Poisson system in a 
convex bounded domain Q with specular reflection on the boundary: 

(1.1) ft + v ■ V x f + V x <j) ■ V v f = , ieOcR 3 ,veR 3 ,t>0 

(1.2) A(j> = [ fdv , x e n , t > 

(1.3) (t, x) = h (x) , xGdn , t>0 

on x 

(1.4) / (0, x, v) = f (x, v) x£tl , veR 3 

(1.5) f (t,x,v) = f (t,x,v*) x£tl , v£R 3 , t>0 

where Q is a convex bounded domain with C 5 boundary, n x denotes the 
outer normal to d£l and 

(1.6) fo(x,v)>0 

Here / (t, x, v) denotes the distribution density of electrons, (j) (t, x) is the 
electric potential. The function h in (jl.3p will be assumed to be positive 
and satisfy the following compatibility condition: 



(1.7) / /o (x,v) dxdv = / h{x)dS x 
Jn Jon 

We will also assume that /o is compactly supported in Q x IR 3 . In (|1.5p 
and in the rest of the paper we use the following notation. Given (x, v) in 
d£l x R 3 we define: 

(1.8) v* = v - 2n x ■ v 

where from now on n x is the outer normal vector to d£l at the point x. 



Key words and phrases. Vlasov-Poisson, Pfaffelmoser method, bounded convex do- 
mains, global existence, specular reflection. 
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The density of particles at a given point x is given by: 

(1.9) p(t,x) = [ f(t,x,v)dv 

In the case of f2 = M 3 , the solutions of the system (|1.1|) - (|1.4|) are globally 
defined in time for general initial data, as it was proved in [13] as well 
as in [11] using different methods. However, in the case of domains with 
boundaries the mathematical theory of well-posedness for the solutions of 
the Vlasov-Poisson system is not so complete as in the case of the whole 
space. It was proved in f6] that classical solutions for the problem (11. lft - 
(I1.5P may not exist in general without the nonnegativity assumption (|1.6j) 
if Q is the half-space M 3 ^. On the other hand, it was also proved in [6] that 
even with the assumption (|1.6j) the derivatives of the solutions of (|l.lj) - (jl.5)l 
cannot be uniformly bounded near the boundary of due to the fact that 
a Lipschitz estimate for the characteristics in terms of the initial data is not 
possible. 

One of the main technical difficulties that must be considered in order to 
solve (jl.ip - (jl.5p . even for short times, is a careful study of the evolution of 
the characteristic curves associated to (jl.ip that remain close during their 
evolution to the so-called singular set, that is defined as follows: 

(1.10) F={(x,v) G n x M 3 : x G dtt, v G T x dtt) 

where T x d£l C M 3 is the tangent plane to dQ at the point x. 

Notice that the projection of such characteristic curves in the domain Q 
bounces repeatedly at the boundary dQ. 

For the case of a half space £1 = Mi the global existence result was shown 
in [7] if the initial data /o is assumed to be constant in a neighbourhood of 
the singular set. The method there is to adapt the high-moment technique 
in |llj . Recently, we have developed in [10] a new proof of global existence 
modifying Pfaffelmoser's idea (cf. |13j). 

In |9j the case of a general convex bounded domain was considered and 
solutions of the linear approximate problem of (|l.ip - (|1.5p were constructed 
under the assumption that the initial data /o is constant near the singular 
set. The absorbing condition was assumed for the distribution density / at 
the boundary to obtain global existence of solution to the full VP system. 
Global existence results in Q = Br (0) were also obtained in the paper, 
under the same assumption, for a class of radially symmetric data that rule 
out possible singular behaviors at the origin. 

Actually, several of the technical difficulties that arise in the study of the 
characteristics near the singular set had been already addressed in [7] in the 
particular case Q = K5_. On the other hand the results in [9] provide tech- 
niques for the problem (jl.ip - (|1.5j) in more general domains. However, the 
assumption of /o being constant near the singular set imposes some restric- 
tions on the initial data, but its main consequence is to make it possible to 
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ignore the evolution of the characteristic curves that are close to the singular 
set. 

The main contribution of this paper is to show how to adapt Pfaffel- 
moser's ideas and to introduce geometric methods to the problem of general 
bounded convex domains with curvatures in order to prove global existence 
in time for the solutions of (|l.l|) - (|1.5j) . It turns out that the effect of the 
geometry of the domains modifies in a stronger way than that of the electric 
field the dynamics of the characteristics. One of the key ideas consists in 
approximating the dynamics of the characteristic curves that are close to 
the singular set by means of the dynamics of a hamiltonian system whose 
trajectories are constrained to the boundary dQ. This approach allows us to 
include easily in the estimates the effects of the curvature of the domain. We 
will then be able to adapt the ideas in [6], [9], [13] to prove global existence. 

The plan of the paper is the following. In Section 2 we introduce a coordi- 
nate system that makes it easier to study the trajectories near the singular 
set and that we will use in the rest of the paper. In this Section we also in- 
troduce a suitable flatness condition for the initial data fo near the singular 
set that will make it possible to obtain solvability for the initial value prob- 
lem for the VP system. In Section 3 we describe an iterative procedure that 
defines a sequence of functions {f n } whose limit as n — > oo yields the desired 
global solution of the problem. In Section 4 we obtain suitable estimates 
for the solutions of the so-called "linear problem" that is the system (jl.ip . 
(jl.4p . (jl.5p with prescribed eft. In Section 5 we prove that the sequence {f n } 
converges to a limit function that is defined as long as a suitable functional 
Q (t) is bounded. Section 6 contains some standard energy estimates for the 
solutions. Section 7 provides a proof of the boundedness of the functional 
Q (t) adapting the ideas of Pfaffelmoser for this problem, in order to deal 
with the geometrical complexity of the domain. This concludes the proof of 
the Theorem. 

2. Preliminary notation and statement of the main result. 

2.1. A more convenient coordinate systems near the singular set. 

By assumption dQ is a C 5 surface and we will parametrized it locally using 
a set of coordinates (^1,^2) • Let us denote as x\\ (/xi,/i 2 ) the point of <9f2 
characterized by the values of the parameters (^1,^2) . We will denote as 
71(^1,^2) the outer normal to <9Q at the point a?n (^1,^2) • 

The Implicit Function Theorem shows that for 5 > sufficiently small we 
can parametrize uniquely the set of points dfl + B$ (0) C M 3 by means of 
the unique values (fj,±, fj,%, x±) solving the equation: 

(2.1) x = xy (/ii,/x 2 ) - a?j_ra(/ii,ju 2 ) 

Given x G dQ + B$ (0) we will represent any vector »eK 3 as: 



(2.2) 



v = uii (hi, fj, 2 ) - v±n (fJ.1,^2) 
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where V[|(/i 1} ^ 2 ) 6 T^^^ (90) , v± £ M and (/Ui,/x 2 ) are as in (|2J])- 
Moreover, we will represent v» = v» (jii, ^2) using the two coordinates 
(^1,^2) defined by means of: 

(2.3) vu = wiu\ + W2U2 

where {^1,^2} are the basis of TL ,( w ,/i 2 ) (90) given by: 

dxii (u-i , no) 

(2.4) n t = L^il^ , i = 1,2. 

The system of coordinates (^1, ^2, x±,wi,W2,v±) will provide a more con- 
venient representation of the set of points in the phase space Hxl 3 that 
are close to the singular set T defined in (jl.lOp . The form that the original 
equation (jl.ip takes in this new set of coordinates is given in the following 
Lemma: 

Lemma 1. The equation 111. 1\) can be rewritten for (x, v) £ [90 + B$ (0)] x 
M 3 , and using the set of coordinates (/ii, ,u 2 , x±, w±, W2, v±) in the form: 

(25) — + V Wl d ^ 1 v — + V<7 — + F— 

dt 1 + kiX 1 9 'Ui dx± 1 dwi dvi 

i=l ' i=l 



where: 



(2.6) «, sE ,-^-± r J^fL, FsE±+ £A 

1 + kiX± ^ 1 + fc 7 x_i_ ^— ' 1 + k 7 x_l 

where kj are the principal curvatures, bj are the coefficients e and g from the 
second fundamental form according to the notation in [15] and Y l - ^ are the 
Christoffel symbols of the surface 90. The vector E = V ' x cj) has been written 
in the form 

(2.7) E = Ekm +E 2 u 2 - E ± n(m,fi2) 
where u\,v,2 are as in pO). 



Proof. The proof of this result is just a standard lengthy change of variables 
that makes use of the classical Gauss- Weingarten equations (cf. [15], page 
124). □ 



Remark 1. The choice of coordinates h2.1\) implies that the coefficient from 
the second fundamental form that is denoted as f in [15] is identically zero. 



Remark 2. Notice that since the domain O is convex, and due to $1.6\) we 
have F < 0. 
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2.2. Compatibility conditions for the initial data. In order to obtain 
classical solutions of (ll.l|) - (ll.5l) we need to impose some compatibility con- 
ditions on the initial data fo (x, v ) at the reflection points of dfl x M 3 (cf. 
[S], [9]). These conditions are the following: 

(2.8) f (x,v) = f (x,v*) 

(2.9) v 1 - [vi7o {x, v*) + V^/o (», v)] + 2E X (0, x) V^/ (x, v) = 



where -E 1 (0, x) is the decomposition of the field E (0, x) given by (|2.7j) 
and V a , V^" are the normal components to d£l of the gradients V x , V„ 
respectively. 



2.3. Flatness condition. The usual way of dealing with the impossibility 
of obtaining smooth solutions for general initial data fo near the singular 
set consists is assuming that fo is constant near such a set (cf. [7] as well 
as [9]). More precisely we will assume that fo S C satisfies the following 
flatness condition near the singular set T: 



(2.10) fo (x, v) = constant , dist ((x, v) , T) < 5o 

for some 5 > small. 

We need to introduce some functional spaces for technical reasons. We 
define for fi £ (0, 1) : 

ll/lloWfixH") = sup il L/m , i 7 Jim +11/1U(QXK3) , V = (V„V t 

V ; (x,v),(x',v')enxR 3 \ \X - X\ +\V -V \ J \ ) 

CY (U x R 3 ) = {/ G C x # (tt xl 3 ): / compactly supported, ||/|| c i,M(fixR3) < °°} 

\Vxf(t,x)- V x f(t',x')\ 

C^Y ([0,T]xfl) = _ SU P r , | r _ r /|M 

+ ll/llc([o,T]xn) + ll/*llc([o,T]xn) 



| V a / (t, x, - V a / x',v)\ + | V„/ (t, x, v) - V v f {t\ x',v')\ 

— _SUp / |M _i_ I /|M ~"~ 

x,i'efi, t,*'e[o,T] |x - x | -+- |« — v \ 

+ ll/llc([0,T]xf2xIR 3 ) + H/*llc([0,T]xQxK 3 ) 

We define the spaces C ([0,T] X O) , C([0,T] xlixR 3 ) as the spaces of 
continuous functions bounded in the uniform norm. 

Remark 3. Theorem\J\ that is the main result of this paper can be proved 
under a more general condition than \2.10\) . namely under the a vanishing 



6 



HYUNG JU HWANG AND JUAN J. L. VELAZQUEZ 



condition similar to the one used in (^\) in the half-line case. More precisely, 

Theorem [JJ is valid under the assumption: 

(2.11) 

<C(x ± +vl) d , 9 

We will explain at the relevant points the modifications that should be needed 
in the proof. However, we have focused mostly in the details of the proof 
under the more stringent assumption \2.10\) for simplicity. 



l/o(x,w)|+ 



1 

x± + v\ 



dfo . . 
-Qx-^ V) 



+- 



1 



1 



X± + V i v± 



dfo f \ 



2.4. The main result: Global existence Theorem. The main result of 
this paper is the following: 

Theorem 1. Let f G C]^ (O xK 3 ),/ >0 for some < n < 1 and 
satisfy i2.10\) and suppose that h G C 2 ' M (dQ) satisfies \1. 7| ) and h (x2, £3) > 
0. Then there exists a unique solution f G ^p\xv) 

((0,00) x n x R 3 ), G 

Ct'x' X ([0> 00 ) x ^) i / or some < A < /i, o/i/te Vlasov-Poisson system 
11. 6\) with compact support in x and v. 

Remark 4. Theorem [JJ cou/c/ 6e derived using similar arguments for the 
case that the function h depends on time and that ^ is smooth enough. 
The arguments would require minor changes but we will just give the details 
of the argument in the case that ^ = for simplicity. 

3. Iterative procedure 

The usual procedure of proving the existence of solutions for Vlasov- 
Poisson models consists in obtaining such a solution as the limit of a sequence 
of functions f n that are defined by means of an iterative procedure. More 
precisely, we define: 

(3.1) f° (t, x, v) = f (x, v) , t > 0, x G n, v G M 3 

(3.2) 

ft + v ■ V x f n + V^- 1 ■ X7 v f n = , xencR 3 ,veR 3 ,t>0 

(3.3) A0 71 - 1 = p n - x (x) = [ f n ~ l dv , xe!l, t>0 

QAU-1 

(3.4) — i = h , x G <90 , f > 

on 

(3.5) / n (0,x,u) =f (x,v) ie!l , «el 3 

(3.6) f n {t,x,v) = f n (t,x,v*) x€dn , v£R 3 , t>0 
for n = 1, 2, .... We assume that /o, /i satisfy (jl.6j) . (|1.7|) as well as (|2.1U|) . 

We will also use the notation: 

(3.7) E n = V0 n 
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Our goal is to show that the sequence f n converges as n — > oo for all 
< t < oo. To this end we need to show as a first step that this sequence is 
globally defined in time for each n > 0. 



4. Linear problem. 

In order to show that the sequence {/"} is well defined we first study 
the well-posedness of the problem (jl.ip . (jl.4p . (jl.5p under the assumption 
that the field E = V</> is given and satisfies suitable smoothness conditions. 
Closely related results have been obtained in [B] for the half-line case where 
geometric complications are not present. 

For further reference we define the evolution of the characteristic curves 
associated to (II. ip . (II. 5p . More precisely, given the field E = V x (j) we define 
for each (x, v) £ fixl 3 the generalized characteristic curve (X (s; t, x, v) , V (s; t, x, v)) 
by the following differential equations: 

dV 

(4.2) — = E = V x ^ 

as 

(4.3) X (t;t,x,v) = x , V (t; t, x, v) = v 

as long as X € fL We extend this definition to arbitrarily long times as- 
suming that at the times s = s* when X n (s*;t,x,v) S dQ the velocity V 
bounces elastically at the boundary, i.e.: 
(4.4) 

V Us*) + ;t,x,v) = lim V (s;t,x,v) = (V ({s*)~ ;t,x,v))* = ( lim V{s;t,x,v) 

s—>s*,s>s* \s—>s*,s<s* 

where (•)* is as in (|1.8p . 

Theorem 2. Assume that E £ C°£f ([0,T] x ft) for some fi £ (0, 1) , and 
E ■ n = h > at dtt. Suppose that f £ C]^ (O x M?) , f > for some 
\i > 0. Then there exists a unique f G Cj.g.'Jj ([0,T] x J2 x M 3 ) , solution to 
the linear Vlasov-Poisson system for some < A < \i. 

Moreover the function f satisfies: 

(4.5) / > 

(4.6) / f(t,x,v)dxdv= / /o (x, v) dxdv , t € [0,T\ 



We will introduce a new coordinate system that will be convenient to 
study the dynamics of the characteristic curves for bouncing trajectories. 
Suppose that (jii, fj, 2 , x±,wi,W2,v±) are as in (|2.ip - (|2.3p . We then define two 
new coordinates (a (t, n\, fj,2, x±, u>i, W2, v±) , ft (t, fix, fi2, x±, wi, W2, v±)) as 



8 



HYUNG JU HWANG AND JUAN J. L. VELAZQUEZ 



follows: 



(4.7) 
(4.8) 

(3 (t, m, H2,X±,Wl,W 2 ,V±) 



a (t,fii,fi2,x±, wi,w 2 ,v±) = 



2ttH (t, hi, fj, 2 ,x±,wi,w 2 ,v±) = it 1 



F (t, fii, H2,0,Wi,W2) x±, 



v± 



2a J 



where the function H will increase by one at each bounce and F is as in 
()2.6p . Since v± changes from —\[2a to \[2a in each bounce, it follows that 
[3 is continuous along characteristics. Notice that (3 is just a coordinate 
that indicates the specific point in the surface {a = constant} where the 
trajectory lies. It does not have the specific meaning of an angle, although 
we have normalized their variation by 2n between bounces. Its functional 
form has been chosen only for convenience. In all the following we will write 
for simplicity F (t, 0) instead of F (t, fJ>i,fJ> 2 , 0, w\, W2) and we will drop the 
dependence of H on the variables H\, fj,2, x±,wi, W2,v± if there is no risk of 
confusion. We can then write: 



o 



x± 



v± 



F(t,0) 



1 



1 



(3 - 2ttH (t) 



IT 



2a 1 



(3 - 2ttH (t) 



7T 



Making the change of variables (t, tii, fj,2, x±, w\, W2, v±) — ► (t, m, H2, a, w\, W2, (3) 
we transform the system ()2.5[) as follows: 

Lemma 2. Under the assumptions of Theorem{l\we have that F (t, fj,±, fj, 2 , 0, wi, W2) < 
—Eq < in d£l in any time interval < t < t* . Moreover, in the coordi- 
nate system (t, fj,±, H2,a,Wi_,W2, f3) the problem \2. 5\) can be reformulated in 
[dtt + B 5 (0)] x R 3 as: 



df Wi df | df 

4-^1 + hx± dm ^ 1 dwi 

i=i i=i 



dt 



+ 



+ 



v x (F (t,xx)-F(t,0))-xx < Y,[j 



Wi dF (t, 0) dF (t, 0) 
+ a, 



TTV 1 



{2a 

irxxvx 

(2a) 3/2 

(4.9) 
= 0, 



3/2 



F(t,0) + 



2vrF {t,0)F(t,xx)x ± 



+ kiXx d/i 



dwi 



(2a 



,3/2 



A / Wi dF (t, 0) OF (t, 0) 

2.1— — + (J < 



. i=l 



+ ^xx dm 



dwi 



d[3 



df_ 

da 
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of the convexity of Q which implies that the term Y^=i i+k x ^ n P-^P is 



Proof. The inequality F (£, //i, A*2j 0,t«i, W2) < — £0 < is a consequence 

negative near dil. On the other hand the continuity of the function h in 
(jl.3p implies that Fj_ = —h < —Eq < whence the result follows. The 
derivation of the equation (|4.9|) follows from (|4.7p - (|4.8p . □ 

Remark 5. We can rewrite lU.Sl) as: 



df | iQj df A 5/ 

1=1 4 = 1 



+ 



+ 




, x . , ^ <9F(i,0) 0F(t,O)' 



ttF (t, 0) + 2vrF (t, 0) [F (t, xj_) - F (t, 0)] x ± 




(2a) 



3/2 



7rx ± w_L I A / Wi dF (t, 0) dF (t, 0) ' 

(4.10) 
= 0, 



9/ 

d[3 



Remark 6. Notice that the dynamics of the tangential part to dQ of the 
characteristics at the singular set is given by the equations: 

dui dwi 

—— = Wi , — — = Oi 

dt dt 

Remark 7. The compatibility conditions A2.8\) , 12. 9\) imply that the original 
data fo written in the variables (t, pL\, fj,2,ct,Wi,W2, (3) is a C 1 '^ function. 
Moreover, the estimate \2.10\) implies that /q is constant for < a < C5q 
for some C > fixed. 

4.1. Velocity Lemma. The following result has been proved in [6], [9] in 
a slightly different manner. 

Lemma 3. Given 5 > fixed we define T$ = ([dQ + Bs{0)] nU) x R 3 . 
Suppose that E satisfies the regularity assumptions in Theorem^ Then, the 
characteristic equations can be solved in the interval t £ [0, T] 

for any (x,v) E Q x R 3 . Moreover, there exist positive constants C\, C2 
depending only on T, fo, ||VF|| loo q t j qi/2^\ such that, for any (x,v) £ 

Tg the following estimate holds: 
(4.11) 

Ci (v 2 ± (0) + x x (0)) < (v 2 ± (t) + x ± (t)) < C 2 (v\ (0) + x± (0)) , < t < T 

Proof. Using the estimate F < —eq < in Lemma [2] as well as (|4.7p , it 
follows that a is equivalent to v\ + x±. Due to the boundedness of $7 it 
is enough to prove this result for small values of v\ (0) + x± (0) , i.e. for 
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points that are close to the singular set. Using (I4.10p we obtain along the 
characteristic curves: 

da , „, NX / Wi 8F(t,0) dF(t,0)' 
~ v± (F (t, x ± ) - F (f, 0))-x± { V —i + <Ti 



dt ~ ~ I 4^ \ 1 + k{X± dm dw 



Notice, however that keeping this term would not change the essence of the 
argument if ^ is smooth as indicated in Remark HI 

Using our regularity assumptions on E we obtain the estimate: 



da 



<C\\VE\\ cl/2(n) \v ± \(x ± ) l/2 + Cx ± 



dt 

where C depends depends only on h and the geometric properties of Oft. 

v 2 

Since -k- + x± < Ka for some positive constant K depending on h, dVl. It 
then follows that: 

(4.12) ^ <Ca 
v ' dt ~ 

Therefore: 

Cia(0) < a{t) < C 2 a(0) , 0<t<T 
for some C\, C2 depending on /o, ||V.E'|| c -i/2(q) , whence (I4.11j) follows. □ 

4.2. Well-posedness of the linear problem. 

Proof of Theorem The proof of this Theorem just follows by integrat- 
ing the equation along the characteristics, combined with the reflecting 
boundary condition. The existence and uniqueness of solutions to (II. ip - 
(jl.5p are not immediate due to the bounces of the characteristics at the 
boundary d£l, but the well-posedness of this evolution has been obtained in 
[9]. We will define the sequence of functions f n by means of the evolution 
of the characteristics associated to (|4.ip - (|4.4p . Then: 

(4.13) / (t, x, v) = f {X (0; t, x, v) , V (0; t, x, v)) 

In order to show that this procedure defines the function / globally in 
time we need to show that for each (x, v) G SI x M 3 the curves defined by 
means of (|4.ip - (|4.4p intersect the boundary dil x M 3 at most a finite number 
of times, and in particular they never intersect with the singular set. This 
fact is a consequence of the Lemma El since for any trajectory starting at 
SI x ]R 3 at time s = t we have a > at time s = t, and therefore a remains 
bounded below and above during the evolution of the trajectory in the inter- 
val s £ [0, t] . Moreover, the characteristics starting in the region {a < C5q} 
where /q is constant, remain in a set of the form {a < C\ (T) 5o} for all 
< s < T, and the characteristics starting in the region {a > C5o} where 
/o is not necessarily constant remain in a set of the form {a > C2 (T) <5o} 
for < s < T. In the first set / =constant. In the second one we have 

that 



is bounded by -S=, where the bound on the number of bounces is 
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uniformly bounded by Finally we notice that, since E £ Cx' 1 ^, classi- 
cal regularity estimates for the solutions of ODEs show that the functions 
X (s; i, x, v) , V (s; t, x, v) are C 1 '^ with respect to the variables (x, v) as long 
as there is no bounces. Moreover, if a trajectory intersects dQ x M 3 our reg- 
ularity assumptions on d£l imply also C 1 '^ regularity with respect to (x, v) 
for the values of the coordinates X where the function X (s; t, x, v) intersects 
9f2, as well as for the time s = s (t, x, v) when such intersection takes place 
(cf. [7j, [9]). Therefore, for trajectories bouncing a finite number of times, 
the functions X (s; t, x, v) , V (s; t, x, v) can be written as the composition of 
a finite number of C 1 '^ functions with respect to the variables (x,v). This 
proves that the function / defined by means of f)4. 13|) is Holder with respect 
to (x, v). The uniqueness of the solution is due to the fact that the solution 
is uniquely detemined by the evolution of the characteristic curves. This 
concludes the proof of the result. □ 

Remark 8. The proof of this Theorem is where there would be a relevant 
difference if the assumption §2.11\) had been used instead of $2.10\) . Indeed, 
if \2.11\) had been assumed, the number of bounces would not be uniformly 
bounded for < t < T. However, it is possible to argue as in [6] to show 
that Theorem [3 holds replacing $2.10\) by §2.11\) . In our setting the way of 
proving this would be to rewrite the characteristic equations associated to 
(ITW as: 

(4.14) 
djM 
dt 

dwi 

If 

da 

lit 

dp 7rF n (t,Q) 2ttF" (t,0) (F n (t,x ± )-F n (t,0))x ± 

— = = 1 3 h hi {m,w i: a, Py/a) 

dt y/2a ( 2a )2 

where the function h^ are smooth in their arguments. The Velocity Lemma 
implies that, for a given trajectory the order of magnitude of a does not 
change in a time interval < t < T. Suppose that we consider trajectories 
with a of order R < 1. Using the change of variables a = Ra, f3 = -^^P we 
would transform ^4-14\ ) in- a similar system of equations for the variables 
C = (fJ'ii w ii 5, p\ with the nonlinearities and their derivatives bounded. 
Classical regularity theory for this system would imply that, for a of or- 



hi (^i,Wi,a,P\fa) 
h 2 (t, m,Wi,a, P\/a) 

v± (F n (t, x ± ) - F n {t, 0)) - x ± h 3 [m,Wi, a, 0^) 



der one 



t with Co = (^Mi,0) Wi.o, Qo, Poj would be bounded, as well as the 
Holder norms evaluated at the points with a of order one. Returning to the 
original variables a, [3 it would then follow that the worse derivative would 



be 



dp 



dao 



that would be bounded as 



In general, the rescaling for each 
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factor a would include a term of order \/R. Due to the decay of /o, V/o 
near the singular set it would then follows that V/ would be bounded like 
a power law near the singular set. Estimates for the Holder norms of Vf 
would then be obtained using the holderianity of fo at distances of order one 
from the singular set, and using the decay of V/ near the singular set, as 

well as the fact that ^f^ 1 '^ 7 ^^m^ — where a is the largest one 

of the corresponding value associated to (x%vi) or (X2V2) ■ 

5. ON THE CONVERGENCE OF THE SEQUENCE {/ n } . 

5.1. Representation formula for the solutions of the Poisson equa- 
tion with Neumann boundary conditions. The following result is stan- 
dard. We just include it here for further references: 

Proposition 1. Given a bounded domain Q in R 3 with a smooth boundary 
d£l, there exists a Green's function G (x,y) for the Laplacian operator with 
Neuman boundary conditions: 

(5.1) A<f> = p(x) , x en 

(5.2) ^ = h , x G dn 

on 

with the compatibility condition 

(5.3) p(x) dx = h (x) dS x 

Jn Jan 

is given by the following representation formula: 

(5.4) <t>{x)= [ G(x,y)p(y)dy- f G (x,y) h (y) dS y 

Jn Jan 

Any other solution of ll5.1\) - [o^3\) is given by \5.4\j up to an additive con- 
stant. The function G (x, y) satisfies the following estimates: 

(5.5) \G(x,y)\ <^=^ , \V x G(x,y)\ <— ^ 

\x-y\ \x-y\ 

C 

\V 2 x G(x,y)\ < -3 , x,yen 

\x - y\ 

where C depends only on the domain fL 
Proof. This result is well known. See for instance [3]. □ 



5.2. The iterative sequence {f n } is globally defined in time. We 

will need the following auxiliary Lemma that states that given / bounded 
in the space C^'A'^ ([0, T] x Cl x M 3 ) , the corresponding field E defined by 
E = V(j), with (f) satisfying (jl.2j) . (II. 3D satisfies the regularity estimates 
required in the Theorems 
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Given a function g : f2 — ► R, we will denote as 

[9] ,X;x 

We define: 



J0,A;a; 



the seminorm: 



sup 



\g (x) -g(y)\ 

\x - y\ 



(5.6) Q(t) = sup {\v\ | {x, v) G supp / (s) , < s < i} . 

We then have the following result 

. Then, the following 



Lemma 4. Suppose that f G C^.L' V ) ([0, T] x $7 x 
estimates hold: 
(5.7) 



(5.8) 



f(t,x,v)dv 



<C(T) 



, x G x [0, T] 



\E(t,x)\<C(T) 
(5.9) 

|F(t,x)| <C{T) 
(5.10) 



O^)([0,T|xnxR3) 



+ 1 , ieOx[fl,T] 



+ 1 , iefflx[o,T] 



\Vp(t,x)\< / \d x f(t,x,v)\dv<C(T)\\f\\ c ^ ([0)T]xnxR3) , *GOx[0,T] 
(5.11) 

|VS(t,x)|+[VS(t,0]o,A;x+[V a ^(t,-)]o,A;«^ C, ( r )ll/llf7^ .an.Tixnx R 3, , xG^x[0,T] 



(5.12) 
(5.13) 



\(H (t,x)\ < C (T) ||/|| c i;M )([0jT]xnxR3) 



\E t (t,x)\<C(T) 



where C (T) > depends on Q (T) and T. 



x G n x [0, T) 
x g n x [0, T] 



Proof. The estimate (15 j7]) is a consequence of the boundedness of the support 
of / as well as the boundedness of the L°° norm of /. The inequality (|5.8p 
follows using a standard regularity theory for the Poisson equation (cf. [3], 
[12]). The estimate (|5.9p is a consequence of the definition of F in (|2.6p 
and our regularity assumptions on d£l and h. The inequality (15. 1QH is just 
a consequence of the regularity properties of / and the boundedness of its 
support. Then (|5.1ip is a consequence of classical regularity theory for the 
Poisson equation. Similarly we can deduce (I5.12P and (15.13p . Thus the proof 
is complete. □ 
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Proposition 2. Let /i,AG (0, 1), satisfying fi > A. Let /q G Cq' m x R 3 ) , 
fo > safe/?/ ^Q73j) . Suppose G C 1 '" (317) , /i > 0. JTien, i/ie sequence 
of functions f n is globally defined for each x G 0, u G M 3 and < t < 
oo. Moreover we have f n G c££*) ([0,T] x O x E 3 ) for any T > and 
ll/ n Hoo = ll/olloo) / Pn(x,t)dx = f f (t,x,v)dxdv. 

Proof. We argue by induction. If n = 1 we use the fact that \ V4>o\ is bounded 
to obtain that f 1 is supported in the region where \v\ < C (1 + t) . Then, p\ 
is bounded by C (1 + 1) . Moreover, due to the regularity of h, we can apply 
Theorem [2] to show that f 1 is bounded in C^h\ ([0,T] x Q. x M 3 ) for any 
T > 0. Applying then Lemma|4]it follows that the norm llE 1 L,i ; i, M , rr , „, 

11 IIU i; X U^M JXS'j 

is bounded. We can then apply Theorem [2] to prove that / 2 is well de- 
fined in C 1,x for < t < oo. Moreover, the support of f 2 would be con- 
tained in the region \v\ < Chi (t) where hi (t) is a continuous increasing 
function in t and using again Theorem [2] it follows that f 2 is bounded 
in C^'h\ ([0,T] x ft x M 3 ) for any T > 0. Iterating the argument we ob- 
tain that the sequence f n is defined as indicated. Finally, using the fact 
that f n just propagates along the characteristics we obtain the conserva- 
tion of || /"Hoc • The conservation of the total mass just follows integrating 
the equation (jl.ip with respect to the variables x, v whence the proposition 
follows. □ 



5.3. The sequence {f n } converges to a solution of the VP system 
if the sequence {Q n } is bounded. We define the following measure for 
the maximal velocities reached for the distribution f n : 

(5.14) Q n {t) = sup{M | (x,v) Gsupp/ n (s), < s < t} . 

Proposition 3. Under the assumptions of Theorem^ suppose that Q n (t) < 
K for n > no, < t < T. Then, f n -> f in C^'\ ([0,T] xSJxI 3 ) as 

n -> oo with < A < jj,, < v < 1 and where f G C^h\ ([0,T]xflx M 3 ) 
is a solution of ( ti.i|) -( f775j) . 

In the Proof of this Proposition we will use some auxiliary Lemmas. The 
first one similar to the one in Theorem 6.2, p 309 of [7]. 

Lemma 5. Suppose the assumptions on TheoremUl are satisfied and that 
Q n (t) <Kforn>n , < t < T. Then: 

(5.15) \E n (x,t)\<C(T) 

(5.16) [E n (-,t)) cl{n) <C(T) , for any 7 G (0,1) 

forn > no+1, < t < T, where C (T) depends only on K, H/o|lL°o(nxR 3 ) ' T. 
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Proof. Estimate (15. 151) follows from classical regularity theory for the Pois- 
son equation as well as the fact that the density p n can be estimated in L°° 
in terms of only Q n (t) and the initial data. 

Indeed, the uniform boundedness of Q n (t) implies that we can estimate 
p n in the interval < t < T uniformly in n. Therefore, regularity theory 
for the Poisson equation implies that E n is bounded in W 1,p (Q) for any 
1 < p < oo uniformly in n. Classical embedding results then imply that E n 
is bounded in C 7 (0) for any 7 £ (0, 1) and the result follows. □ 

We now prove the following basic Lemma that shows that the boundedness 
of Q n implies the boundedness of / in the norm C^.'hs ([0,T] x Q x R 3 ) . 

Lemma 6. Suppose that Q n (t) < K for n > no, < t < T. Then: 
(5.17) linic^ao.Tlx^^Cm , 0<t<T 

for A < fj,, n > no + 1, where C (T) depends only on K and T. 

Proof. Indeed, the estimates (j5. 15[) . ()5. 16j) imply that E n is uniformly bounded 
on the Holder spaces Cl for < 7 < 1. Notice that choosing 7 > \ we can 
derive the estimates in the velocity Lemma (cf. [3|) uniformly in n. In par- 
ticular this implies that a > C5$ uniformly in n for t S [0, T] . 
We can write the characteristic equations for (|4.10p as: 

(5.18) 

= h (t,m,Wi,a,(3) 

~^ = v± {F n (t, x±) — F n (t, 0)) - x ± f 3 fawi, a, (3) 

dp 7rF n (t,0) 27rF"(t,0) (F n (t,x ± )-F n (t,0))x ± 

-77 = 7^= 1 3 \- H(fH,v>i,a,P) 

dt V2a (2a) 5 

where the functions fi, fy, f\ depend only on the regularity properties of 
dQ, and the boundary value h, and therefore are smooth. The function fi 
depends also on the field E. 

We can now take the Holder derivative of (|5.18p with respect to the initial 
data Co = {pi.OiWifl^ao, Pq) . Let us write also Q = (fJLi t ,Wi t ,a, P) and: 

r 1 \g(x ,v ) -g(x' ,v' )\ 



|(a;o,«o)-(a:| ) ,« / )|<l ICo ~ Col 



Since, a is uniformly bounded below and using also the uniform bound- 
edness of Q n it follows that for any A < u [C]a ^ satisfies an inequality of 
the form: 

d 



dt 



<c[C], 



■■Co 
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Therefore [C]a-c — ^ where C (T) depends only on K and T. Us- 
ing the formula /q (xq,vq) = f n (t,x,v) , we obtain a uniform estimate for 
f n (t, ■, •) in , (fi x M 3 ) for any < A < fi. It then follows that p n is uni- 
formly bounded in C% (fi) and using classical regularity theory for the Pois- 
son equation it then follows that E n (i, •) is uniformly bounded in Ci' A (fi) . 
Using again the lower bound for a we can now apply standard regularity 



results for ODEs to (|5.18p to obtain that 



9Co 



is uniformly bounded for 

A;Co 



< t < T. Therefore f n is uniformly bounded in C}^ v -, (fi x M 3 ) . Using then 
the equation (|3.2p we obtain uniform estimates for /" in Cfc, x (fi x M 3 ) . 
This implies (I5.17j) and the Theorem follows. □ 

Remark 9. The result in LemmalM can be proved if $2.1Q\) is replaced by 
t2.11\) using the rescaling argument in Remark^ 

The following simplectic property is a standard consequence of the fact 
that the evolution of the characteristics curves is hamiltonian. 

Lemma 7. Moreover, let us denote as (X (s; t, x, v) , V (s; t, x, v)) the so- 
lution of the characteristic equations \4.1\ )- \4.!fy . For any s G [0, T] the 
transformation 

(x, t) -> (X (s; t, x,v) ,V (s; t, x, v)) 
is simplectic. In particular: 

(5.19) dX (s; t, x, v) dV (s; t, x, v) = dxdv 
We finally conclude the Proof of Proposition [3j 

Proof of Proposition We recall that the sequence f n has been defined 
by the iteration (|0) - (j577]) . 

We claim that the sequence is a Cauchy sequence in L 1 ([0, T] x fi x M 3 ) . 
To prove the claim, let f n+1 and f n be consecutive elements of the sequence 

{f n Y 

(5.20) + v ■ V x f n+1 + V x r ■ V v f n+1 = 0, 

(5.21) f? + v V x f n + V^"" 1 • V v f n = 0. 

Substr acting (I5.21|) from (|5,20p yields 
(5.22) 

(r +1 - p) t +v-v x (r +i - r)+v s ^-v„ (r +1 - r) = (v,^- 1 - v^o-v^r. 
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By integrating f|5.22 j> along the trajectory (X (a) , V (a)) with X (t) = x and 
V (t) = v we get: 

(5.23) 

(r +1 - r) (x, v ,t) = { r +i | t=0 - ru) (* w , no)) 

+ / (V^™- 1 - V x 0") (a, X (a)) • V„r (a, * (*) , V (a)) da 
(V^ 1 - V x n ) (a, X (a)) • V„/" (a, X (a) , V (a)) da 



where: 

Applying the representation formula (|5.4j) and the estimates (|5.5p to com- 
pute the difference (V x (/) n ~ 1 — V x (j) n ) (s,x) we obtain: 

,,_ , n _i _ , nW \\^n [ \P n (v) ~ P"' 1 Ml j 
\(V x (f> -V x (p ){s,x)\<C J ■ -J L dy 

Jsi | a; — y\ 

Integrating (15. 23ft over the phase space (x, v) we obtain, applying Fubini's 
theorem: 



f a+1 {t) f n (*)|| L l( nxM 3) 
ft 



< 




J JQxR :i 



[V x <t> n - X - Vx<t> n ) (s, X(s))\ |V„/ n (a, X (a) , V (a))| dvdxda 



t 

n-l 



<C/ / X"(y,a)| /0 "(a,y)- /0 ™- i (s,y)|^, 
JO isi 

where: 

^0^) = // TV77 -|V„r(s,^(s),^( S ))|dX(a)^(a) 

i JqxR 3 \X (s) - y\ 

- — ^—i 1 \V v f n (s,x,v)\dxdv. 
'sixR 3 \x — y\ 

In the last identity we used the Liouville principle (|5.19p . We now estimate 
K n (x,t) as follows: 

# (x, t) = / ■ -5 -dy + / ■ - 2 v -dy 

J\y-x\<r \x-y\ J\y-x\>r \x — y\ 

<o||v./-ft»)li rW) + ' V,r y" IiW) . 

We choose r to optimize the right hand side of the above inequality, namely: 

l|V v / n (i)lliia^ 
r||V w /»(<)llL ? ( i i, = - ^F^- 
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where from now on we use by shortedness the notation Lf {L\) = L°° (J^L 1 (R 3 )) , L\ (Lj) 
L 1 (O x R 3 ) . Thus we have 

K n (t,x) < C \\V v f n (t)\\ 2 £(TM H V -/ n ^ " V3 



Since f n is bounded in W '°° and the supports are bounded uniformly in 
n (due to the global bound on Q n (t) , for < t < T) it is easy to see 
that ||V„/ n an d ll^o/ 71 (*) 11^00(1/1) are uniformly bounded in n 

and hence -fT n (i, x) is uniformly bounded. 

Thus we obtain 
(5.24) 



\\f n+ Ht)-f n (t)\\ L1{nxm <C(T)J*\\P(s)-f 



n, ~ 1 fslll 



1 ( s )|L 1 (filxK 3 ) ds 



where C = C (T) depends only on T, Q (T) , and the initial data. Notice 
that (15. 24ft implies by iteration that: 

(5-25) l|/" +1 (*)-/ n (i)L 1(f ,x R 3)<^ 

for some 6 < 1, and < i < £q, if eo is sufficiently small depending only on 
T. Then (|BT24l implies: 

|| r +l {t) _ f n ( t )\\ < C (T) df+C (T) /" I)/"' (,) - / 

Therefore we obtain that (|5.25p is valid for t £ [eo,2eo] changing C\ 
if needed. Iterating the argument, it then follows that {/ n } is a Cauchy 
sequence in the space L°° ([0, T) , L 1 (O x M 3 )) . 

Once we know that f n is a Cauchy sequence in L 1 ([0, T] X f2 X M 3 ), 
we can show that the sequence is Cauchy in C^'}^\ ([0, T] X X R 3 ) for 
any < A < /i, < < 1 arguing by interpolation. Indeed, using 
(|5.17p with A replaced by A satisfying A < A < fi and interpolating be- 
tween L 1 ([0, T] X O X R 3 ) and TV 1 ' 00 ([0,T] X f2 X R 3 ) we obtain that / n is 
Cauchy in W l,p ([0, T] x $7 x M 3 ) for any p > 1. Using Sobolev's embed- 
dings we can obtain that / n is a Cauchy sequence in C v ([0,T] x x M 3 ) 
for any < v < v. Interpolating then in Schauder spaces between C 1,A 
and C u we obtain the desired convergence. In order to check that / 6 
Cji'Jj ([0,T] x O x M 3 ) , we use (EH) to obtain: 

/ n (*) = h- f[v v,r ( S ) + v^™- 1 ( s ) • v„r ( s )] ^ 

JO 

Passing to the limit in this equation as n — ► oo it follows that: 

/ (*) = fo ~ ! [v ■ V./ (s) + V x <p (s) ■ V v f (s)] ds 
Jo 

and this implies the desired differentiability for /, and the Proposition fol- 
lows. □ 
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5.4. Convergence of the sequence Q n to Q. 

Proposition 4. Let Q n , Q be as in |5.i^[ ), i t 5. 6\) respectively. Suppose that 
max {sup n>no Q n (t) , Q (i)} < A' /or < t < T. Let us assume also that 
f n -> f in C%h\ ([0,r] X n X M 3 ) /or any < A < fjt, < v < 1. T/ien 
limn^oo Q n (i) = Q (i) uniformly on [0, T] . 

Proof. The characteristics starting in a (0) > C5q remain during their evo- 
lution in the set {a (t) > C8q\ due to Lemma [3j Therefore, these charac- 
teristics remain separated from the singular set and we can estimate their 
difference as n — > oo as it was made in [10] . Indeed, for these characteristics 
the number of bounces is uniformly bounded in n in the interval < t < T. 
Moreover, the times where these bounces take place for the functions f n 
converge to the corresponding times for the bouncing times for the charac- 
teristics associated to / and using the fact that E n — > E as n — ► oo we obtain 
also that the corresponding characteristic curves associated to f n converge 
to the ones associated to / between bounces, and due to the boundedness of 
the number of such bounces, it follows that the functions V n (s, t; x, v) con- 
verge uniformly to V (s, t; x, v) as n — > oo. Since Q n (t) is the maximum value 
of \V n \ associated to characteristics (X n (s,t; x,v) , V n (s, t; x, v)) which at 
s = lie in the support of /o, it follows that Q n (t) — > Q (t) , and the result 
follows. □ 



5.5. Prolongability of uniform estimates for the functions f n . 

Proposition 5. Suppose that for some T > there exist K > and uq > 
such that for any n > uq and < t < T we have Q n (t) < K. Then, there 
exists eo = Eo (K, H/olloo) > such that for < t < T + eo and n > no the 
following estimate holds: 

Q n (t) < 2K 



Proof. Notice that: 
(5.26) \p r > 



f n dv 



<[|/0[loo 



Using the representation formula (15. 4|) for the solutions of the Poisson 
equation in f2 with Neumann boundary conditions in Proposition [1] we ob- 
tain: 

|v<n < c H/olL (Q n (t)) 3 + J Pn dx + \\h\lcx,, 

On the other hand, using (|4.2j) and (I5.14|) we obtain the following inequal- 
ity for t >T: 

Q n+1 (t) < Q n (T) + C H/olL £ (Q n (s)f ds + C(t-T) 
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where C > is just a numerical constant independent of n and Q n . Using 
the assumption on Q n we obtain: 

Q n+1 (t)< i^ + C H/olL f (Q n {s)fds + C(t-T) , t>T 

Jt 

Defining R n (t) = max {Q^ (t) : uq < i < n} , we obtain: 

i?"+i(t)<ir + C||/ || oo f* {R n+1 (s)) 3 ds + C(t-T) , t>T 

Jt 

Let us select en = „r„. , ¥~ „■» n ■ It then follows, using a Gronwall type of 
C[8||/ |Ui<r 3 +i] 

argument that: 

Q n (t) < 2K , n > n , < t < T + e . 

□ 



6. Energy estimate and consequences. 

The following energy estimates are standard for the Vlasov-Poisson sys- 
tem (cf. [5]). We state them here for further reference. 

Proposition 6. Suppose that f is a solution of hl.l\) - [Lb^) on the time 
interval t E [0, T] . Then: 

(6.1) 4- ( I v 2 fdxdv + I (Ef dx] = 

at \JnxR 3 Jn J 

Proof. It is similar to the proof of the analogous result in the whole space 
(see [5], page 120). The only difference is that we need to take into account 
the contribution of some boundary terms. More precisely, using (jl.ip we 
obtain, after some integration by parts: 

dt \JnxR 3 J JnxR 3 JdQxR 3 

(6.3) 

f E 2 dx)=-2[ V x( f>-vfdxdv + 2 [ dS x 
dt \JnxR 3 J JnxM. 3 Jan \ an xJ t 

+ 2 / (j>fv ■ n x dS x dv 
JdUxR 3 

The boundary term in (|6.2p vanishes, since Jqq xR 3 v 2 v -n x fdS x dv = for 
each x E dfl, due to the specular boundary condition (jl.6p . On the other 



hand, the two last terms in (|6.3p vanish due to the fact that = h is 
independent of t and L 3 fv ■ n x dx = for each x E dVL. Adding then (|6.2p . 
(EH) we obtain (Ell. □ 
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Proposition 7. Suppose that f is a solution of $1. 1\ )- (T75\) defined in < 
t < T with f(0,x,v) = fo(x,v), where /o is as in fl.6\) . There exists C 
depending only on T and on the regularity norms assumed for fo in Theorem 
[7] such that: 



(6.4) 
(6.5) 



sup / v f (x, t) dvdx < C 
o<t<T Jn 



sup 

0<t<T 



|p(V)ll 



< c 



(6-6) ll/(*)lliP(n X R3) = H/olLp(nxR3) , for all 1 < p < oo. 

Proof. This result is just a consequence of ()6.1|) and its proof is standard in 
kinetic theory. See for instance [5]. Finally (|6.6p follows multiplying (jl.ip 
by (/) n ~ j and integrating by parts using the specular boundary conditions. 

□ 



7. PFAFFELMOSER'S ARGUMENT: GLOBAL BOUND FOR Q (t) . 

In this section we show that the function Q (t) can be bounded in any 
time interval < t < T and therefore that the corresponding solutions of 
(jl.ip - (|1.5j) can be extended to arbitrarily long intervals. The estimate of 
Q (t) will be derived using the ideas of Pfaffelmoser (cf. [13]) in the case 
of bounded domains with purely reflected boundary conditions at d£l. The 

main content of the result is a uniform estimate for Q (t) as long as / is 
defined. 

Arguing as in the derivation of (|5.26p we obtain the following estimate: 

(7-1) IMI^H/ILQW 3 

where p is in (|1.2[) . 

The main result of this paper, that we will prove using Pfaffelmoser 
method, is as follows: 

Theorem 3. Let f G C 1 ^ (£1 x M 3 ) with < (J, < 1. Suppose that f £ 
C t t X v) ([0, T] x f2 x M 3 ) is a solution of (Ejp-(E1|) with A G (0, 1) , < T < 
oo. There exists a (T) < oo depending only on T, Q (0) and ll/ollci.^nxiR 3 ) 
such that: 

(7.2) Q(t)<a{T) , < t < T . 

7.1. Bounds for Q (t). Suppose that (^X (s) , V (s^j is any fixed character- 
istic curve such that: 

(1(0), V(0)) Gsupp/o. 

The basic idea in the Pfaffelmoser's method consists in deriving estimates 
for E = X7(f). Using the estimates in Proposition Q] it follows that: 
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(7.3) 



ds 



t-A 



E[s,X(s) 



< C 



c 



t 

t-A 
t 

t-A 



ds 



ds 



f(s,y,w) 



flxl 



y-x(s) 

f{t,x,v 



r dydw + C \\h\\ A 



T dxdv + C 



X(s)-X(s) 

where we are using the following change of variables: 

(X (s) , V (s)) = (y, w) - (x, v) = (X (t) , V (t)) 

We also note the measure preserving property (dydw = dxdv) that is due 
to the fact that the evolution of the characteristics is hamiltonian away from 
the boundary and that the measure dxdv is also preserved by reflection on 
the boundary. 

Pfaffelmoser's method is based on the idea of splitting the region of in- 
tegration in (|7.3p into three different sets that are usually termed the good, 
the bad, and the ugly. Fix Q > 0, that will be precised later. In the rest of 
the argument we will use two numbers A and P defined by 

(7.4) 



(7.5) 



A = 



5 > small 
cpP 

Q 2 



where cq is small, but fixed number (independent on Q). 

Given Q, A, P we define the good, the bad, and the ugly respectively by 



G = { (s, y, w) e [t - A, t] x n x R 3 : \w\ < P, w-V (t) 



w-V + (t) 



<P 



B = {(s,y,w) E [t- A,t] x n 



w - V (t) 



w - V + (t) 



U = {{s,y,w) G [t- A,t] x O 



w - V (t) 
where 



> P, 



w - V + (t) 



i J : \y-X(s)\ <e , \w\ >P, 

>_ P] 

R 3 : \y-X(s)\ >e , \w\ > P, 

>_r} 

R 



1 



v-V(t) 



if the characteristic curve X (s) does not intersect d£l on the interval s G 
[t — A, t] , and otherwise 



i? 



1 



1 





- + - 




v — V (t) 




v-V+ (t) 
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Here 
where 



V + (t)=V* (so), 



s = sup e [t - A,t] : X (s) € dtt\ . 

The change due to the field can be estimated as: 

Lemma 8. Under the assumptions in Theorem we have the following 
estimate: 

ft 



(7.6) 



\E(X(s),s)\ds < CA (Q(t)) 4/3 + l 



t-A 



, A<t<T 



Proof. It is just an adaptation of a similar estimate in [5] , page 122. The last 
constant term in (|7.6p is a consequence of the Neumann boundary condition, 
i.e., E = h at the boundary (cf. (|1.3p ). □ 



We now estimate the change of V (s) due to geometry of the domain and a 
basic geometric property that relates the change of V (s) with the curvature 
will be investigated: The idea is as follows: 

(i) The total length that the characteristic moves, including reflections is 
bounded as CQA. 

(ii) The change of the normal vector in a distance of order li is bounded 
by Ci % . 

(iii) Therefore, the change of the angle of the vectors V (s) , V* (s) with 

respect to the previous reflection is bounded by (cQli + J^ +1 \E\ ds^ . (The 

change is the difference of outcoming vector with respect to the next out- 
coming one, and incoming vectors with respect to the incoming ones). 

(iv) The total change of these vectors is then bounded by CQ 2 A + 
f*_ A \E\ds < C (Q 2 + Q 4 / 3 ) A < CQ 2 A, where the constant C depends 
only on the maximum of the curvature. 

We now give its explicit formulation and the rigorous proof. To this end 
we first recall basic ingredients such as Formulas of Frenet in differential 
geometry (see [15] p. 18, p. 94 for reference) that we will combine with 
Lemma [2] in this paper. 

Lemma 9. Let T and N be the unit tangential and normal vector on d£l 
respectively. Let k be the curvature along a curve and be the normal 
curvature in the direction dx of the line of curvature. Then we have 

dN + Kjydx = 0, dT = nNds, 

where s is the arc length. 



Proof. See [15] p. 18, p. 94. 



□ 
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In view of Lemma [2] the equations describing the evolution of the charac- 
teristics in a geometrical form are: 



(7.7) 
(7.8) 

(7.9) 
(7.10) 



dt 

dwi 
~dt 

dx± 
~dT 
dv± 
~df 



W; 



1 + kiX± 

v±Wiki 



E 



1 + kiX± 
v± 

3=1 



2 rywfwj 

■ , i 1 + k i x± 



+ kjX± 



Note that these equations work only near the boundary dtt, but inside 
the domain the maximum displacement of the characteristics is QA that 
can be made small. 



Lemma 10. Let (X (s) ,V (s)) and IX (s) ,V (s) ) be two characteristics 
and let Q be as in \5. 6\) . Then we have 



mm 



{|y(s)-y(t)|,|y(s)-F + (t)|}<CAQ 2 (t) , s£[t-A,t] 



mm{\V(s)-V(t)\,\V(s)-V + (t)\} <CAQ 2 (t) , s e [t - A,t) 

Proof. Estimates equally apply to V and V and so we only give a proof for 
V. Unlike the whole space case, we have to take care of the possible sign 
change of v± at the bounces , i.e. if it becomes zero coming from the region 
v± < 0. Using the equations (|7.8|) and (|7.10j) . it follows that 

<CQ 2 , a€[t-A,t] 

and this implies 

(7.11) \ Wi (s)-Wi(t)\ <CQ 2 (t-s), 



d(\v±\) 
ds 



(7.12) \\v±(s)\-\v ± (t)\\ <CQ 2 (t-s). 

(Notice that the equation above is valid even if we cross the bounces, and 
that \v±\ is differentiable) . If the number of jumps of the normal velocities 
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is even, then by (17.llh - fl7.12n and by Lemma O we have 
\V(s)-V(t)\ = 



W i U i + V ± N ^J ( S ) ~ WiUi + V±N ^j W 

2 

(m (a) - wi (t)) (t) + («i (s) - v± (t)) iV (t) 

i=i 

2 

+ ki («) I k (*) - «i («) I + |"l («)| 1^ (*) " ^ 



i=l 

< CQ 2 (t - a) . 

In a similar manner, we deduce that if the number of jumps of the normal 
velocities is odd, then 

\V(s) -V + (t)\ < CQ 2 (t-s). 

Thus we conclude the assertion of the lemma. In general, we can split 
the time interval [s, t] into several time sub-intervals in such a way that 
particles are governed completely by the equation (|2.5p near the boundary 
or by the usual Vlasov equation (II. ip away from the boundary on each sub- 
interval. Then the estimates above combined with those in the whole space 
and noticing that the effect Q (t) 4//3 of the electric field alone (Lemma [8|) is 
negligible to that Q (t) of the geometry yield the lemma. □ 



Lemma 11. In the sets B and U, we have 



w 



— < \v < 2\w\ , 



w - V (t) 




w-V + (t) 




-+- 





< 



v-V(t) 


+ 


v-V + (t) 


< 2 


w - V (t) 


+2 


w - V + (t) 



> p, 



Proof. In the sets B and U, we have \w\ > P, w — V (t) 

P. By Lemma ITU] we have either \w — V (t)\ < CAQ 2 or 
CAQ 2 . U\w-V (t)\ < CAQ 2 , we have 



w-V+(t) > 
w-V+(t)\ < 



W-V(t) — | V (t) — XI7| < 



v — V(t) 



< 



w-V + (t) -|V(*)-«;|< v-V + {t) < w-V + {t) 
Since, for small cq 



w-V{t) +|V(t)-«>| ) 
+ \V(t)-w\ 



CQ 2 A <-< 



CQ 2 A < — < 



w-V(t) 



w-V+ (t) 
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we have 



w — V (t) 



w-V+ (t) 



< 



< 



v — V (t) < 2 w - V (t) 



v-V* (t) 



< 2 



w 



In the other case, we similarly obtain 
w - V + (t) 



< 



v — V (t) < 2 w - V + (t) 



\w-V(t) 



< 



v-V* {t) < 2 w — V (t) 



Therefore, we deduce our lemma. 



□ 



We denote = fi\Ui + ^2^2, X± = x± and Vj| = w\Ui + W2U2, V± = v±. 

Lemma 12. // a trajectory (X, V) has more than one bounce in the interval 
[t — A, t] , then we have, for all s £ [t — A, t] , 

\V±(s)\ <CQ 2 (t)A. 

Proof. If a trajectory (X, V) has more than one bounce, then we have 
V± (s) = 0, for some s G [t — A, t] . Since 

d\V ± \ 



(7.13) 
the lemma follows. 



ds 



<CQ Z (t), 



□ 



Lemma 13. Let (X (s) , V (s)) and (^X (s) , V (s)^ be trajectories over [t — A,t]. 
Suppose that 



X± (s ) - X ± (so, 



mm 

se[t-A,t] 



X±(s)-X ± (s) , so G (t- A,t). 



Then either both X± (s ) > 0, X± (s ) > or both X± (s ) = X± (s ) = 0. 

Proof. We prove the lemma by contradiction. Suppose X± (sq) > and 

X± (so) = 0. The case X± (so) > and X± (so) = can be studied in a 

2 



X± (s) ~ X ± (s) 



is differ- 



symmetric way. Notice that the function A (s) 

entiable in the interval (t — A, t) except at a finite set of points. Therefore, 
at the point s = so where the minimum of A is achieved we have: 



X ± (so-) - X ± (so-) 



< 



dX . . d 
ds ds 

where from now on / (s -) = lim s ^ S0;S<S0 / (s) , f (s +) = lim s ^ S0)S>S0 / (s) . We 
then have: 



X± (s ) (V± (so-) - V ± (so-)) < 0, 
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which implies 

(7.14) V ± (s -) < V± (s -) < 0. 

The fact that V± (so— ) ^ is a consequence of Lemma El 
Notice that, since X± (so) = there is a reflection of V± at s = so. On 
the other hand, since Xj_ (so) > , V± is continuous at s = so- Therefore: 

VI (ao+) = VI (s -) , V±(s +) = -V ± (s -) 
Thus we have 

X ± (s)-X ± (s)| = X ± (s)-X ± (s) =X L (s )+(vi(flo-) + Vl(«o-)) (s-s )+o( S -s ) 
as s — > sq, s > sq. Due to (|7.14p . we have X± (s) — X± (s) < X± (sq) = 



*L («o) " X ± (s ) 



for s — sq > sufficiently small, but this contradicts the 



fact that 



X± (s) - X± (s) 



reaches its minimum at s 



result follows. 

We show the following crucial separation property. 



so- Therefore the 
□ 



Lemma 14. (Separation property) In the ugly set U, there exist sq,s\ S 
[t — A, t] such that the following separation holds: 
(7.15) 



X(s)-X(s) 



> C ( Sq + min 



{\v-V(t) 



•sol 



V + (t) 



Sl 



}) , se[t-A,t] 



where C is a universal constant depending only on the curvature of d£l. 

Proof. We separate into two cases: 

Case 1: Both trajectories (X (s) ,V (s)) , (x(s),V(s) \ have at most 

one bounce in the time interval [t — A, i]. Let t — A < t\ < t<i < t be 
possible two bouncing times with t%,t2 corresponding to (X (s) , V (s)) , 

(^X (s) , V (s)^ respectively. We split the time interval [t — A, t] into a max- 
imum of three sub-intervals, namely [t — A, t\] U [ti, t2\ U fo, t] . In the most 
general case, some of these intervals could be empty. Let us describe the 
argument in the most general case, since for a smaller number of reflections 
the argument required is just a minor simplification of it. Pick so and si 
such that 



mm 

se[t-A,ti]u[t 2 ,t] 

min 

se[ti,t 2 ) 



X 
X 



X) (s) 
X) (s) 



X -X 
X -X 



(so) 
(si) 



In the interval [ti,i2] there is no bounces along the trajectories. Then, we 
can argue exactly as in the case without boundaries (cf. [5], pages 128-129) 
to show that: 



X(s)-X(s. 



> C 



V + (t) |s-si|, for s G [ti,t 2 \. 
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On the other hand, the portion of the trajectories X (s) ,X (s) for s G 
[t2,t] might be reflected with respect to the boundary d£l. Suppose for the 
moment that the boundary dfl is flat as in the half space case. The trajecto- 
ries over [t2,t] obtained by reflections together with the original trajectories 
X (s) , X (s) for s G [t — A, t\] yield portion of new trajectories without 
bounces, and satisfying an equation of the form: 

dX 

ds 
dV 

ds 



V 



E 



where J t _ A E (X (s) , s 



ds < CA 



(Q(t)) 4/3 + l 

actly as in the case of the whole space (cf. [5]), to estimate the difference 
between the trajectories, and since the reflection with respect to the plane 
x\ = is an isometry, we finally obtain: 
(7.16) 



. We can argue then ex- 



X(s)-X(s) 



> C 



V(t)-V(t) 



sol 



for sG [t- A,ti] U [t 2 ,t] . 



Now if the boundary d£l is not flat then the change of the normal vectors 
between two reflections can be bounded by CAQ and the corresponding 
change of the vectors V and V can be bounded by CAQ 2 , which is smaller 
than cqP for sufficiently small cq. Therefore the change of these vectors is 

small compared to V (t) — V (t) in the ugly set and the inequality (|7.16p 

above is valid for the case of domain with curvature. 

Case 2: At least one of the trajectories has more than one bounce in 

[t — A, t] . Let (^X, V^j have more than one bounce in [t — A, t] . 

We first consider the case 

1 



(7.17) 



\wi (s) 



Wi 



> 



i.e., the tangential part of V 
whole interval [t — A, t] . Let 

min luj (s) 
se[t-A,t] 



V (s) - V (s) , for all s G [t - A, t] , 
- V dominates along the trajectory in the 



fa 0)1 = I Mi ( s o) - fa (so) I 



Note that the tangential part of the trajectory, X\\ is C . Consider the 
equation (|7.7|) and (|7.8|) for the tangential components of the position and 
the velocity: 



dwi 
~dT 



Ei 



v±Wjki 
1 + kiX± 



^ 1 + kjx± 

j,£=l 3 ± 



o(q 4 / 3 )+o(q 2 ), 



djM 
ds 



W; 



1 + kiX± 



Wi 



(s ) + O (AQ 2 ) + O (xl H) = Wi (s ) + O (AQ 2 ) , 
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where the second term comes from the change on the velocity Wi due to the 
field and to the geometry and constants depend only on the geometry of the 
domain. We also used the fact that 

\x±\\wi\ =0(AQ 2 ). 

Then, integrating the tangential part fa — fa of the difference of (X, V) 
and ( X, V ) , we get: 



(m - fa) (s) = (m - fa) (s )+(wi (s ) - m (s )) (s - s )+O (AQ 2 ) (s - s )) 
By Lemma [10] and fjT. 17f) , we have 



\wi (s ) - Wi (s )| > - V (s ) - V (s ) 



> 



1 . r 

— mm < 

2 X 



v-V(t) 



V + (t)\}+0 (Q 2 (t)A). 



Since | (/i« — /tj) (s)| 2 is a C 1 function it folllows that at the point s G 
[t — A, i] where |(/x, — /tj) (s)| 2 attains the minimum we have 

( s_s °)^ — Ai) («)| 2 ) 

= [/Xi (s ) - fli (s )] • K («o) - Wi (s )] (s - s ) > 0. 

Note that the inequality above takes into account the possibility that the 
minimun of \(fii — fli) (s)| 2 can be achieved at the end points sq = t — A,t. 
Since 



O {Q 2 it) A) = c P, 



we deduce 



10** — Ai) ( s )\ > ^minj 



v — V(t) 



v-V 



(t)\} 



s - s ) 



for sufficiently small c®. Using the definition of the ugly set we have |(/i, — fa) (s)\ > 
e . Thus (I7T5D follows. 

We now consider the complementary case, i.e., there is s £ [t — A,t] such 
that 

1 



> 



V(s)-V(s) 



V± (s) - V ± (s) 
By ()7.13p . we have 

|Kl(s)|-|Vi(s)|| > \\V± (I)|-|Vl(I)| I -CQ 2 (t)A. 
On the other hand, note that 



\V±(s)\ 



V±(s] 



V ± (s)-sgn(V ± (s)) V ± (s) 
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Using Lemma \12\ it follows that 



V ± (s)-sgn(V ± (s)) V ± (s) 

V± (a) - V± (s) + V± (s) - sgn (V ± (s)) V± (s) 



< 
< 



V\ (s) - V\ (s) 



+ 2 



V±(s) 



V±(s)-V ± (s) +CQ 2 (t)A. 



Similarly, we get 
V\ (s 



sgn(V ± (s)) V ± (s) 
V± (I) - V± (I) + VI (I) - sgn (V ± (s)) V± (s) 



> 
> 



V\ (s) - V\ (s) 



V± (s) 



V±(s)-V ± (s) -CQ 2 (t)A. 



Thus, we obtain, for all s £ [t — A, t] 



(7.18) V ± (s) - V ± (s) 



> 



V\ (s) - V\ (s) 



CQ 2 (t) A 



> 



V(s)-V(s) -CQ 2 {t). 



> - min{ v-V(t) , v-V + (t)}- CQ 2 (t) A 

> - min{ v-V(t) , v - V + (t) } - c P 

1 - - P 

> - min{ v-V(t) , v - V + (t) } > — . 

Using Lemma \12\ it follows, choosing cq in (|7,5p sufficiently small, that: 

P 



(7.19) 



\Vl{s)\> 



for all s £ [t — A, t] . Taking into account (|7.6p it follows that V±(s) changes 
sign, by reflection, at most once in the interval s £ [t — A, t] if cq is suf- 
ficiently small. Combining Lemma [T2l (|7.18p . an d (|7.19p it follows that 
V±(s) — V± (s) changes sign at most once for s € [t — A, t] . Indeed, V±(s) 
is small compared to \V± (s)| > in the interval [t — A,t], and V± (s) 
changes sign only once at most. Suppose that V±(s) changes sign at s = so- 

Since X± (s) > 0, it follows that sign (V± (s)) = sign (v± (s) - V± (s) 

sign (s — sq) , for s G [t — A, t] . We have 



X± (s) - X ± (s) 



so 



V± (r) - V± (r) 



dr. 
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Then, using (|7.18j) . we have: 

X ± (s)-X ± (s)\ >imin{ Vlh- Y(l) . !'(/)- V + (l) }|.s-.s„| 



Since in the ugly set £/, 
is complete. 



X± (s) - X ± (a) 



> £q, we obtain (|7.15|) . The proof 

□ 



Proof of Theorem 0. The key point in the proof is to estimate the right- 
hand side of (|7.3|) . Let us assume without loss of generality that Q > 1, 
since otherwise the leading contribution would be CA in (|7.3p . In order to 
make this estimate we separate the contributions of the sets G, B and U : 

/ (y, w, s) 



ds 

t-A JnxIR 3 



y-x(s) 



jdydw 



dsdydw 



G 



...] dsdydw + 

B JU 



dsdydw 



In order to estimate the contribution of the good set we define: 



Po(y,s) = / f(y,w,s)dw 
Jg 

Standard estimates yield 

IIpgIL < cp* 

where from now on C depends on ||/o|loo > no * on ^ £ o 5 R, Q, A. 
Arguing as in the derivation of (|T.6[) we obtain: 



(7.20) 



dsdydw < AP 4/3 



In order to estimate the contribution of the bad set, notice that Lemma [TT] 
implies 

8R I 1 1 



w 





- + - 




w-V(t) 




w-V+ (t) 



Then 



(7.21 

C 



f(y,w,s) 



ds 



t-A 



y-x(s) 



R 



dsdydw < C 



ds 



t-A 



\w\<C 



\w\<Q \W\ 



where by assumption 



V(t) 



Eodw 



dw < CRA log (Q) 



> 1, since otherwise the corresponding charac- 





- + - 




w — V (t) 




w-V+ (t) 



teristic does not have effect in the variation of Q. 
Lastly we estimate the integral over the ugly set: 



:!2 
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The aim is to estimate 

f{y,w,s) 



u 



dsdydw 



f(x,v,t) 



y-X(s) J u X {s)-X{s) 

Using (|7.15p we can estimate this integral as 

f(x,v,t) 



rdsdxdv 



C 



u 



Eq + min - 



v-V(t) 



\s - sol , 



v-V+ 



(t)\) 



rdsdxdv 



s — Sl 



The integration with respect to time s can be estimated as 

I* ds 

Jt-A 



(eq + m i n | v — V (t) \s — sq\ , v — V + (t) \s — s i|j) 



< 



< 



C 
£o 

c_ 



1 



1 



v-V(t) 


- + - 


v-V* (t) 


1 


- + - 


1 


v-V{t) 


v-V* (t) 



Then 



r CI' C 

(7.22) J [...} dsdydw <— J v 2 f (x,v,t)dvdx < —A 

Combining the estimates (j?T20|) . (f7T2Tj) . and ([T722]) for the Good, Bad, and 
Ugly set, we obtain 

J^Je(s,X (s)) \ds<CA (V 3 + iilog (Q) + 



C A \PV 3 + R log (Q) + 



Q 4/3 ' 
RP 



Choosing R = Q 1 ' 6 , P = Q 3 / 4 " 5 we obtain 



t-A 



b a(»: 



/3 



where /3 < 1. 
Therefore 



Q(t)-Q(t- A) 
A 



< CAQ 



<C(Q(t)f, 



and a standard iteration yields Q (t) bounded in any interval < t < T. 
Thus this completes the proof. □ 
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Proof of Theorem [7J In order to prove the existence of a solution / of 
(|l.ip - (|1.6j) globally denned in t we will show that the sequence of functions 
f n defined by (13. 1 () - (|3 . Tf) converges as n — > oo to a solution of f)3.2f) - (13.6() 
for arbitrary values of t. To this end, it suffices to show that the functions 
Q n (t) are uniformity bounded in each compact set of t G [0, oo) .The desired 
limit property would then follow from Proposition [3j 

To this end, we define L (t) = sup n Q n (t) . We have that L (•) is increasing 
in t. We denote as T max the time where: 

lim Lit] = oo 

f- * max 

Our goal is to show that T max = oo. Let us assume that T max < oo. We define 
£o = £o (2<r (r max ) , H/olloo) , where the function £q (') is as in Proposition 
[5] and the function a (T) is as in Theorem [3j Notice that the functions 
Q n (i) are uniformly bounded for t G [0, T max — by definition of T max . 

Therefore, Proposition [5] implies that / n — > / in C^'/' » for < i < T max — 
4j, < v < 1. We can use Proposition [4] to prove that Q n (t) — > Q (t) 
for t G [0, T max - si] . In particular, lim^oo Q n (i) = Q (i) < cj (T max ) , for 
t = T max — 4^. Therefore Q ra (t) < 2a (T max ) for n > no with no large. Then, 
Proposition [5] implies that the sequence Q n (t) is uniformly bounded for 
t G [0,T max + 4J 1 ] , whence L (t) is bounded as t — > T max . This contradicts 
the definition of T max and concludes the proof of the existence of a solution 
of (|l.ll) - (ll.6p in Cj'i v \ for some < A < fi and for < t < T as asserted in 
Theorem [TJ In order to prove uniqueness we argue as in the proof of (|5.24l) 
to obtain, that two C^'/j | ^ solutions of (|l.ip - (|1.6p with the same initial and 
boundary satisfy 

ll/i (t) - h {t)\\ L x < C (T) f \\h (s) - h (8)\\ L i ds, 

Jo 

Therefore f\ = fa. This completes the proof of the Theorem. □ 
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